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ABSTRACT 

Although delayed detonation models of thermonuclear explosions of white dwarfs seem promising for re- 
producing Type la supernovae, the transition of the flame propagation mode from subsonic deflagration to 
supersonic detonation remains hypothetical. A potential instant for this transition to occur is the onset of the 
distributed burning regime, i.e. the moment when turbulence first affects the internal flame structure. Some 
studies of the burning microphysics indicate that a deflagration-to-detonation transition may be possible here, 
provided the turbulent intensities are strong enough. Consequently, the magnitude of turbulent velocity fluctu- 
ations generated by the deflagration flame is analyzed at the onset of the distributed burning regime in several 
three-dimensional simulations of deflagrations in thermonuclear supernovae. It is shown that the correspond- 
ing probability density functions fall off towards high turbulent velocity fluctuations much more slowly than a 
Gaussian distribution. Thus, values claimed to be necessary for triggering a detonation are likely to be found 
in sufficiently large patches of the flame. Although the microphysical evolution of the burning is not followed 
and a successful deflagration-to-detonation transition cannot be guaranteed from simulations presented here, 
the results still indicate that such events may be possible in Type la supernova explosions. 
Subject headings: Stars: supernovae: general - Hydrodynamics - Instabilities - Turbulence - Methods: nu- 
merical 



1. INTRODUCTION 

The question of whether deflagration-to-detonation transi- 
tions (DDTs) may occur in type la supernova (SN la) explo- 
sions has been a puzzle since the scen ario of delayed d etona- 
tions was suggested in the first place (lKhokhlovlll991") . If a 
supersonic detonation burning mode occurs at all in SNe la, 
such a transition is inevitable because a prompt detonation 
fails to prod uce the observed intermediate mass elements 
(lArnettlll969h . For these to be synthesized, burning must par- 
tially take place at low fuel densities. A pre-expansion of the 
white dwarf (WD) material before burning is only possible 
for subsonic propagation of the thermonuclear flame. Conse- 
quently, the standard picture is that of an accreting WD which 
getting close to the Chandrasekhar limit becomes unstable to 
a thermonuclear runaway and ignites a subsonic deflagration 
flame. In contrast to the shock-driven detonation, flame prop- 
agation is mediated by microphysical transport in this mode. 
However, this scenario is intricate because the energy released 
in burning may expand the WD so rapidly, that the slow flame 
is not capable of burning sufficient amounts of material be- 
fore it has dropped to densities at which burning cannot be 
sustained any longer This way, the explosion would be far 
too weak to be consistent with observations. Therefore, the 
flame needs to accelerate and two possible mechanisms for 
this have been suggested. 

The first is well-founded on known physical principles and 
is based on the fact that the buoyant rise of the burning 
bubbles generates strong turbulent motions due to Rayleigh- 
Taylor and Kelvin-Helmholtz instabilities. The flame inter- 
acts with eddies of the resulting turbulent cascade and, by 
means of surface area increase, accelerates significantly. This 
turbulent deflagration scenario has been rec ently studied in 
detail i n three-dimensional simulations (e.g. 'R einecke et alJ 
l2002bt iGamezo et alJ llOOl lR5pke"&H illebrandll l2005h . It 
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has been shown that it may lead to explosions of the WD star 
meeting the gross properties of observed SNe la ( Ropke et alJ 
,2007a) . However, some issues remain unsolved in this 
model. First, in its current implementations it can explain 
only the weaker events. Second, it may in some configu- 
rations leave behind unburnt material in the ce ntral parts of 
the ej ecta, which disagrees with observations (Kozm a et alj 
I2OO5I) . Third, it has problems explaining the high- velocity in- 
termediate mass elements seen in the spectra of SNe la. 

These weaknesses of the models may possibly be cured if a 
second way of flame acceleration, a DDT, takes place in later 
stages of the explosion. One-dimensional parameterizations 
of this "delayed detonation model" showed best agreement 
with observations assuming this transition to occur once the 
fuel density ahead of the flame has dropped to ^lO^gcm"^. 
Unfortunately, such a transition is hypothetical as of yet, since 
no convincing mechanism has been identified that would work 
r obustly in a SN la. 

iNiemeyer & Woosleyl (Il997h pointed out that the only fun- 
damental change in the flame properties is found at the tran- 
sition from the flamelet regime of turbulent combustion to the 
distributed regime. In the first, which holds for most parts of 
the explosion process, the interaction of the flame with tur- 
bulence is purely kinematic. The flame front as a whole is 
corrugated by turbulent eddies, but its internal structure re- 
mains unaffected. At a first glance this seems surprising, 
since Reynolds numbers Re{L) at scales of L = 10^ cm typi- 
cal for the situation in SNe la are of the order of 10'"^. The 
corresponding Kolmogorov scale, 77 L[Re{L)]~^l'^, down to 
which the turbulent cascade extends, is ^ 10""* cm — much 
smaller than the flame thickness. However, assuming Kol- 
mogorov scaling, the velocity fluctuations (denoted as v' in 
the following) decrease with the cubic root of the length scale 
under consideration within the inertial range of the turbu- 
lent cascade. This constitutes the so-called Gibson scale, at 
which the laminar flame speed is comparable to v'. Below the 
Gibson scale, the flame burns through turbulent eddies faster 
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TABLE 1 

Parameters of the three-dimensional deflagration SN Ia 
simulations used for the analysis. 



name 


spatial 


resolution 


central 




coverage 


(computational 


density 






grid cells per 


[lO'gcm-^] 






dimension) 


Simulation I 


full star 


1024 


2.9 


Simulation II 


full stai' 


640 


2.0 


Simulation III 


octant 


256 


2.9 



than they can deform it, and therefore they do not affect the 
flame shape considerably. However, since the laminar burn- 
ing velocity and the flame thickness depend on the fuel den- 
sity O^immes & Wooslev 1992), the Gibson scale decreases in 
the course of the explosion while the flame structure becomes 
wider. Therefore, turbulent eddies will eventually penetrate 
the internal flame structure and mix heated material, fuel, and 
ashes. 

It turns out that this onset of the distributed burning 
regime takes place at fuel de nsities favored for DDTs i n 
one-dimensional simulations jNiemever & WooslevI 119971) . 
Studying the micro physics of burning in this regime, 
iLisewski et al.l (l2000t) concluded that triggering a detonation 
may be possible for sufficiently strong turbulence. The veloc- 
ity fluctuations should be close to 10^ cms"' (corresponding 
to ^20% sound speed), but at this time the available two- 
dimensional supernova simulations seemed to exclude such 
high turbulent velocities. A re-analysis of this scenario is un- 
derway (Wooslev 2007) and arrives at similar conclusions. If 
a mechanism providing a deflagration-to-detonation transition 
driven by the flame itself existed, the parameter deciding on 
triggering a detonation would be the strength of the mixing of 
the flame structure, which is determined by the magnitude of 
v'. This value can be estimated from three-dimensional nu- 
merical simulations of the deflagration stage of the explosion 
and is the focus of the present paper. 

2. APPROACH 

To determine v', data from several three-dimensional simu- 
lations of the deflagration stage in thermonuclear supernovae 
are analyzed. These simulations apply th e methods desc ribed 
in detai l by ReineckeJtiD (fT999i l2002al) . iRopkd (l2005h . and 
ISchmidt et al.i (i2006i) . The key features include flame track- 
ing via a level-set approach, a turbulent subgrid-scale model, 
and a moving computational grid to follow the expansion of 
the WD. 

Simulation I dRopke et alj|2007ah was carried out on 1024^ 
grid cells and comprised the full star, as did simulation II 
which was set up on 640^ grid cells. Only an octant of the 
WD was accounted for in Simulation III with a 256^ cells grid 
assuming mirror symmetry with the other octants. The ex- 
ploding WD had an initial central density of 2.9 x lO^gcm""* 
in simulations I and III, while for simulation II a central den- 
sity of 2.0 X lO^gcm"^ was chosen. These parameters of the 
setup are summarized in Table [T] All initial WD configura- 
tions were composed of a mixture of equal parts, by mass, 
of carbon and oxygen. The flame was ignited in multiple 
kernels spherically distributed around the center of the star 
an d partially overlapp ing, similar to the scenarios presented 
bv lRopkeetan(l20()6l) . 

The quantity under consideration here is v' experienced by 
the flame in the distributed burning regime. Of course, this is 
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Fig. 1 . — Fits to the histogram of v' in simulation I at t = 0.80 s for 1 < 
p[10^gcm"^] < 3 (top two panels), and the probabihty of finding v' larger 
than a given value according to eq. (2) for different times in the same density 
interval (bottom panel). 

not directly resolvable in simulations carried out on the scales 
of the WD star Therefore, an estimate of the turbulent veloci- 
ties at the flame is obtained from the subgrid-scale turbulence 
model. It provides the value of v' at the scale of the compu- 
tational grid. Since the size of the grid cells is a dynamical 
quantity in the implementation, it is reseated to a length scale 
of 10'' cm for comparison. To this end, a Kolmogorov spec- 
trum was assumed, which need not necessarily apply to the 
actual scaling of turbulence at the cons idered spatial range ; 
however, in the simulations presented bv lZingale et alj ( l2005h 
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TABLE 2 

MAXIMUM TURBULENT VELOCITIES AT A SCALE OF 10* cm, FIT PARAMETERS ACCORDING TO EQ. (Q}, PROBABILITY OF FINDING v' > 10* cms"' 
(P(10*)), ESTIMATED FLAME AREA Aest, SIZE OF THE PATCH OF THE FLAME WHERE v' > 10* cms"' (Aesl/'(10'*)), AND NUMBER OF COMPUTATIONAL 
CELLS TAKEN INTO ACCOUNT IN THE ANALYSIS FOR SIMULATION I IN THE DENSITY RANGE 1 < p[10''gcm"^] < 3. 



f Is] 


v' no'' 


cm) 


-ao [10-^] 


ai [10-'] 


-612 


P(10*) 




Acst [cm-] 


AestPdO*) 


cells 




[10^ cms 
















[cm2] 




0.70 


6.70 




0.356691 ±0.0004 


7.12569 ±0.0006 


12.4931 ±0.007 


5.062E 


-07 


7.85£ + 16 


3.98£ + 10 


55,341 


0.75 


8.38 




2.97325 ±0.006 


5.98318 ±0.001 


11. 1650 ±0.02 


1.702£ 


-06 


5.35£ + 17 


9.11£ + 11 


274,015 


0.80 


11.5 




15.3360±0.010 


5. 12493 ±0.0004 


9.90905 ±0.007 


2.179£ 


-06 


1.50£ + 18 


3.28£ + 12 


574,769 


0.85 


10.1 




8.92170 ±0.02 


5.44507 ±0.001 


10.3121 ±0.02 


5.016£ 


-07 


2.52£ + 18 


1.27£ + 12 


736,593 


0.90 


7.05 




8.92137 ±0.006 


5.55688 ±0.0004 


9.47029 ± 0.007 


8.966£ 


-09 


2.98£ + 18 


2.67£ + 10 


614,195 


0.95 


6.20 




5.18888±0.01 


5.87899 ±0.001 


9.71264 ±0.02 


1.691£' 


-09 


2.89£ + 18 


4.89£ + 09 


443,823 


1.00 


4.24 




2.20126 ±0.005 


6.39985 ±0.001 


9.98804 ±0.02 


6.35 IE 


-11 


2.28£; + 18 


1.45£ + 08 


263,740 



TABLE 3 

Selected parameters for different density intervals and 

SIMULATIONS. THE NOTATION IS THE SAME AS IN TABLe[2] 



f [s] 


v' ( 1 0* 


cm) AestP(10*) [cm-] 


cells 




[10' cms 


-'] 




Simulation 1, 1 < p[10'gcm-'] < 2 


0.70 


5.37 


6.72£ + 05 


5723 


0.75 


7.47 


1.05£+11 


80,438 


0.80 


11.5 


1.66£ + 12 


276,683 


0.85 


10.1 


1.20£ + 12 


443,841 


0.90 


6.63 


1.82£ + 10 


423,654 


0.95 


6.20 


9.74£+09 


355,975 


1.00 


4.24 


2.58£+08 


243,480 


Simulation 1, 2 < p[\0' gcra^] < 3 


0.70 


6.70 


6.00£+10 


49,618 


0.75 


8.38 


4.40£ + ll 


193,577 


0.80 


9.30 


9.83£+ll 


298,086 


0.85 


7.79 


3.06£ + ll 


292,752 


0.90 


7.05 


1.78£ + 09 


190,541 


0.95 


4.81 


6.96£+08 


87,848 


1.00 


3.22 


4.92£+07 


20,260 


Simulation 11, 1 < p[10' gcmT^] < 3 


0.75 


5.27 


4.52£+10 


8312 


0.80 


6.65 


1.95£ + 11 


45,902 


0.85 


8.36 


8.37£ + 10 


121,638 


0.90 


8.46 


1.37£+11 


200,764 


0.95 


7.02 


1.80£ + 11 


239,816 


1.00 


6.65 


6.70£ + 10 


230,245 


1.05 


5.79 


2.91£ + 11 


187,809 


1.10 


5.18 


3.38£ + ll 


139,548 


1.15 


4.45 


1.67£ + 10 


90,103 


Simulation 111, 1 < p[10'gcm-3] < 3 


0.75 


5.18 


5.11£+12 


9865 


0.80 


5.67 


6.17£ + 10 


15,073 


0.85 


6.51 


1.39£+12 


18,209 


0.90 


5.75 


3.80£ + 12 


17,166 


0.95 


4.04 


2.61£+11 


12,589 


1.00 


3.00 


2.10£+08 


8098 



and iRopke et alj ( l2007al) . Kolmogorov turbulence is indeed 
recovered for buoyancy-unstable flames. In any case, the po- 
tential error introduced by this procedure is expected to be 
small, since the grid scales at times examined here will not 
deviate significantly from 10^ cm. Another simplification is 
introduced by reducing the data sets by a factor of 2 in each 
direction, leaving out every second cell of the computational 
grid. Although this deteriorates the statistics, it simplifies 
the analysis of the large data sets from the highly resolved 



simulations. The distributed burning regime is expected to 
be reached for fuel densities in the range 1 . . . 3 x lO^gcm"-* 
(IWooslevll2007l: iNiemever & Kerstein|[l997l) . Therefore, we 
measure v' determined by the subgrid-scale model in compu- 
tational grid cells cut by the flame front in this density range. 
This introduces two uncertainties. First, these cells contain 
mixed states of fuel and ashes. Consequently, the values we 
obtain are not equal (but, due to the small density jump over 
the flame, similar to) the values in pure fuel. Although it may 
seem compelling to restrict the analysis to cells sufficiently 
far away from the flame so that they contain pure fuel, this 
would require a distance of at least three cells, and turbulence 
will be considerably weaker that far ahead of the flame. Val- 
ues derived this way would therefore not well represent the 
y' experienced by the flame front. Second, as we are looking 
for maxima of v', the subgrid-scale model will provide only a 
mean value on the size of the grid cells. Of course, there will 
be a maximum of this quantity, which will be determined in 
the following, but due to the intermittent nature of turbulence, 
rare extrema exceeding these values are possible on smaller 
scales in the real situation. 

3. RESULTS AND IMPLICATIONS 

The following discussion focuses on simulation I (the cor- 
responding values are given in Tables|2]and[3]i. Results from 
the other simulations are included in Table [3] and corrobo- 
rate the generality of these results. The possibility of trigger- 
ing a detonation at a c ertain turbulence strength depe nds on 
the density of the fuel jLisewski et al.ll2000t [moslevl l2007). 
To account for this density dependence we consider v' at the 
flame front in two separate density ranges in simulation I: 
(1 < P7 < 2) lO^gcm"* and (2 < < 3) lO^gcm'^ 

The instantaneous maxima of v' in the snapshots of the sim- 
ulations analyzed here fall around 10^ cms"' for all three se- 
tups (cf. Tables |2] and O. They show a similar temporal evo- 
lution passing through a peak value. The maxima of v' are 
slightly lower for less resolved simulations. Since the turbu- 
lent subgrid-scale model should compensate for decreasing 
resolutions, this is not likely due to lower turbulence predic- 
tion. Instead, a lower resolution implies a less frequent real- 
ization of high turbulence intensities on grid. 

3.1. Probability of high turbulent velocities 

The information gained from determining the maximum 
turbulent velocities in the simulation is limited since only 
stochastic realizations are recovered. The robust physical 
quantity that can be extracted from the simulationss is the 
probability density function (pdf) of v'. This pdf can be ap- 
proximated with the normalized histogram of v' determined 
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Fig. 2. — Simulation laXt = 0.80s after ignition. The blue opaque, transparent, or wiremesli surfaces correspond to tlie flame front (implemented via a level-set 
approach). Volumes of high turbulent velocity fluctuations are rendered in red/orange. In the lower regions, the green arrow indicates the location of the maximum 
value of v' found in the simulation. For better visibility, white areas correspond to ash regions and fuel regions are shown in black in a plane intersecting with the 
maximum v' -value in the lower right region. 



as described above. An example is shown in Figure (Tlfto/^j. 

For determining the possibility of detonations, the high- 
velocity tail is the relevant part of the histogram. An estimate 
of the corresponding part of the pdf is obtained by fitting the 
histogram with the exponential of a geometric Ansatz: 



P(v') = exp[flo(v'r'+fl2] 



(1) 



The fit parameters for simulation 1 at different times are given 
in Table |2] Integrating this expression from v to oo yields the 
probability of finding v' greater than a given value v. 



P{v')dv' = 



V exp(a2)r (l/ai,-aov" 
ai(-ao)''''"v 



(2) 



where r(., .) denotes the upper incomplete Gamma function. 
The results of this fitting procedure for an exemplary data set 
are shown in Figure [T] Obviously, Gaussian or lognormal fits 
fail to reproduce the high-velocity trend of the histogram (see 
Fig.[T] middle). 

Figure[T]('fooftomj illustrates the time evolution of the prob- 
ability of finding a turbulent velocity larger than a given value, 
according to equation At earlier times, the probability of 
finding v' > 5 x 10^ cms"' is relatively large, but it decreases 
steeply for higher threshold values. Later on, this decrease 
flattens out somewhat, but the overall probability becomes 
lower. 
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To discuss this behavior, we introduce the probability 
of finding v' > lO^cms"', P(IO^), because according to 
iLisewski et all teoOOh . this is approximately the threshold for 
triggering a detonation. Values of ^(10*^) are given in Table|2l 
For simulation I they follow the temporal trend of v[^^^ and 
decrease steeply towards late times. This evolution might be 
expected. Turbulence decreases in the late stages of the su- 
pernova explosion because the energy injection from buoy- 
ancy instabilities becomes lower and due to expansion (t he ex- 
plosion ultimately establishes homologous expansion; Ropke 
12005). Thus, the probability of finding high turbulent veloc- 
ities on any patch of the flame naturally decreases with time. 
But here we consider only the part of the flame for which the 
fuel densities are in a certain range. The corresponding flame 
area increases first and then decreases, since the fuel densities 
steadily decrease in the explosion process. This area can be 
estimated from the number of computational cells cut by the 
flame front and their sizes, providing the value Aest given in 
Tables |2] and [3] This temporal evolution of Aest makes real- 
izations of a certain high v' more likely at a particular instant 
given by the convolution of the area effect with decreasing 
turbulence strengths. A quantification is provided in Tables |2] 
and [3] with the size of the flame surface Aest^'(10'^) at which 
v' > 10^ cms"'. Obviously, this area needs to exceed a cer- 
tain size in order to be relevant for a DDT. As seen from Ta- 
ble|2l the evolution of Aest^'(lO^) follows the maxima of v' de- 
termined at different times in simulation I. However, for the 
other simulations the trend is less clear. In simulation II the 
density range of interest is reached later due to lower initial 
central density. This delays the generation of turbulence and 
thus burning, energy release, and the expansion of the WD. 
The overall lower values of Aest^'(10'^) derived from simula- 
tion II may be taken as an indication that lower central densi- 
ties decrease the chances of a DDT, but this hypothesis needs 
to be tested in a larger sample of simulations. 

3.2. Location of high turbulent velocities 

In order to establish a self-sustained detonation wave, a 
minimum mass of material needs to be burned in the initi- 
ation process (Niemever & Wooslevl ll997l; iDursi & TimmesI 
12006; Ropke et al. 2007b). Therefore, it is necessary that the 
high turbulent velocities are concentrated in larger patches 
and do not distribute all over the flame in tiny regions. Al- 
though Aest^'( 10*^) exceeds lO'^cm^ in simulations I and III, 
there is no guarantee that a connected spatial structure with 
dimensions larger than 10km exists. Moreover, the particu- 
lar location at the flame where the high turbulent velocities 
are fo und may be critical for the detonation ignition (IWoosleyl 
l2007h . Do they occur at leading or trailing edges of burn- 
ing bubbles, or in small channels of fuel in between large ash 
structures? Such geometrical information cannot be provided 
by the pdf. 

To find answers to these questions we consider the situation 
at 0.8 s after deflagration ignition in simulation I. 

Values of v' exceeding 10*^ cms"' are directly found in the 
simulation (see Table|2]i. The fact that the grid resolution here 
is ^10* cm indicates that such high turbulent velocities are lo- 
calized in finite patches. If the area of high turbulent velocity 
fluctuations determined from the pdf was distributed in many 
very small patches, it would not be directly visible in the sim- 
ulation, due to numerical smoothing. 

Figure|2]shows the location of the regions of high turbulent 
intensity at the deflagration flame. The upper left region of 
Figure |2] demonstrates that flame regions with high turbulent 



intensities are present mainly at trailing edges of bubble-like 
structures. The location of all regions with turbulent velocities 
larger than 10^ cms"' is shown again in the upper right region, 
where the flame has been made transparent in order to make 
these regions better visible. At the instant shown in Figure |2l 
the peak value of v' in the analysis of the simulation is found. 
Its location is illustrated in the lower regions. 

3.3. Possibility of triggering a detonation 

According to iLisewski et al.l (l2000l) . values of v' > 
10** cms"', corresponding to about 1/5 to 1/4 of the sound 
speed at the considered densities, are necessary for a DDT 
(with optimistic assumptions on the DDT mechanism). The 
exact values depend on the fuel density and composition. For 
a fuel composition of equal parts, by mass, of carbon and 
oxygen (this case applies to the simulations analyzed here), 
a threshold of v' > 5 x 10^ cms"' is given for a fuel density 
of 2.3 X lO^gcm"^. For lower fuel densities, this thresh- 
old becomes larger, and at 8 x 10''gcm"-' it is given with 
v' > 8 X 10^ cms"'. The maxima found in the simulations (in 
particular, in simulation I; see Tables |2] and [3]l exceed these 
thresholds. 

The estimates for the total size of the flame patches where 
v' > 10^ cms"' and the fact that high values of v' are found 
at scales resolved in the simulation indicate that the corre- 
sponding regions may be sufficiently large to trigger a deto- 
nation (compare with the va lues of detonator sizes given by 
iNiemeyer & Woo slev' ("1997^ and'R opke et"an (l2007ta) ). 

As noted by [Lisewski et al. (20Q3), flie chances of trig- 
gering a detonation increase with higher carbon fractions 
in the fuel and decrease with higher oxygen fractions. 
iRopke & Hillebrandtl (|2004|) noted that the energy release in 
burning does not depend strongly on the carbon mass fraction 
of the fuel as long as it proceeds to nuclear statistical equi- 
librium. This, however, does not apply to the situation con- 
sidered here. The low fuel densities imply incomplete burn- 
ing to intermediate-mass elements. Therefore, in addition to 
the effect of altering the thresholds of turbulence strength re- 
quired for triggering a detonation, the carbon mass fraction 
may also affect the turbulence strength in the distributed burn- 
ing regime. Additional studies are required to settle this ques- 
tion. 

4. CONCLUSIONS 

The intensity of turbulence was analyzed at the onset of 
the distributed burning regime in three-dimensional SN la 
simulations. Considerably larger values for the maximum 
turbulent veloc ity fluctuations than previously anticipated 
(ILisewski et al.l l2000) were found for all simulations consid- 
ered and the conclusion that these values are typical for three- 
dimensional simulations featuring a state-of-the-art treatment 
of turbulence on unresolved scales seems compelling. The 
histogram of v' features a pronounced high-velocity tail. This 
part is not well fit by a Gaussian pdf but rather by an expo- 
nential of a geometric Ansatz. 

According to our result s, a deflagratio n -to-de tonation tran- 
sition as anticipated by Lisewski et alj (1200 0) is not ruled 
out. But, lacking a microphysical model, they cann ot provide 
certain ty either, and the concerns pointed out by iNiemeyeJ 
d 19991) persist. However, the results presented here encourage 
detailed studies of the microphysics of the distributed burning 
regime at hi gh turbulent velocities, such as those presented by 
ILisewski et al. (2000) and Woosley (2007|) . 

If a DDT occurs in a SN la explosion, then a location at the 
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outer parts of the deflagration flame, but at the trailing edge of 
a bubble-like feature, seems most likely. The rapid decline of 
the pdf towards high turbulence intensities indicates that such 
a transition may be a rare event, possibly realized only once 
or a few times in a supernova. This, however, as well as the 
question of whether a DDT occurs in every event, depends 
on the details of the formation of the detonation and warrants 
further study. 

It may thus be possible that detonations do not form at all 
leading features of deflagration flame but only in a few such 
locations. As a detonati on front cannot cross even tiny regions 
filled with nuclear ash (iMaier & Niemeye3l2006h . it may not 
reach all patches of unburnt material embedded in the com- 
plex deflagration structure. However, in numerical simula- 
tions of dela yed detonations in WD with parametrized DDTs, 
iRopke & Nie mever (2007) found that even a detonation front 
formed in a single spot burns most of the remaining fuel by 



reaching the dense center of the star. Such large-scale su- 
pernova simulations with parametrized DDTs help to address 
the question of whether the outcom e of delayed de tonations is 
consistent with observational data (.Mazzali et al.ll2007) . 



Stan Woosley greatly supported this work with stimulat- 
ing discussions. The simulations have been performed on fa- 
cilities of the Max Planck Society, Garching, Germany, the 
HPCx, Edinburgh, UK (as part of the DEISA project), and at 
the National Center for Computational Sciences at Oak Ridge 
National Laboratory, which is supported by the Office of Sci- 
ence of the U.S. Department of Energy, under contract DE- 
AC05-00OR22725. The work was supported by the SciDAC 
Program of the DOE (DE-FC02-01ER41 176) and the NASA 
theory program (NNG 05-GG08G). 



REFERENCES 



Ai-nett, W. D. 1969, Ap&SS, 5, 180 

Dursi, L. J., & Timmes, R X. 2006, ApJ, 641, 1071 

Gamezo, V. N., Khokhlov, A. M., Oran, E. S., Chtchelkanova, A. Y, & 

Rosenberg, R. O. 2003, Science, 299, 77 
Khokhlov, A. M. 1991, A&A, 245, 114 

Kozma, C, Fransson, C, Hillebrandt, W., Travaglio, C, Sollerman, J., 

Reinecke, M., Ropke, R K., & Spyromilio, J. 2005, A&A, 437, 983 
Lisewski, A. M., Hillebrandt, W., & Woosley, S. E. 2000, ApJ, 538, 831 
Maier, A., & Niemeyer, J. C. 2006, A&A, 451, 207 

Mazzali, R A., Ropke, R K., Benetti, S., & Hillebrandt, W. 2007, Science, 
315, 825 

Niemeyer, J. C. 1999, ApJ, 523, L57 

Niemeyer, J. C, & Kerstein, A. R. 1997, New Astronomy, 2, 239 
Niemeyer, J. C, & Woosley, S. E. 1997, ApJ, 475, 740 
Reinecke, M., Hillebrandt, W., & Niemeyer, J. C. 2002a, A&A, 386, 936 
— . 2002b, A&A, 391, 1167 

Reinecke, M., Hillebrandt, W., Niemeyer, J. C, Klein, R., & Grobl, A. 1999, 
A&A, 347, 724 



Ropke, R K. 2005, A&A, 432, 969 

Ropke, R K., & Hillebrandt, W. 2004, A&A, 420, LI 

— . 2005, A&A, 431,635 

Ropke, R K., Hillebrandt, W., Niemeyer, J. C, & Woosley, S. E. 2006, A&A, 
448, 1 

Ropke, R K., Hillebrandt, W., Schmidt, W., Niemeyer, J. C, Blinnikov, S. 1., 

& Mazzali, R A. 2007a, ApJ, 668, 1132 
Ropke, R K., & Niemeyer, J. C. 2007, A&A, 464, 683 
Ropke, R K., Woosley, S. E., & Hillebrandt, W. 2007b, ApJ, 660, 1344 
Schmidt, W., Niemeyer, J. C, Hillebrandt, W., & Ropke, R K. 2006, A&A, 

450, 283 

Timmes, R X., & Woosley, S. E. 1992, ApJ, 396 
Woosley, S. E. 2007, ApJ, 668, 1109 

Zingale, M., Woosley, S. E., Rendleman, C. A., Day, M. S., & Bell, J. B. 
2005, ApJ, 632, 1021 



